Complete transfer of populations from a single state to a pre-selected superposition of 
states using Piecewise Adiabatic Passage: Experiment. 
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We demonstrate a method of adiabatic population transfer from a single quantum state into 
a coherent superposition of states. The transfer is executed with femtosecond pulses, spectrally 
shaped in simple and intuitive manner, which does not require iterative feedback-controlled loops. 
In contrast to non-adiabatic methods of excitation, our approach is not sensitive to the exact value 
of laser intensity. We show that the population transfer is complete, and analyze the possibility of 
controlling the relative phases and amplitudes of the excited eigenstates. We discuss the limitations 
of the proposed control methods due to the dynamic level shifts and suggest ways of reducing their 
influence. 
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INTRODUCTION 

Population transfer from one energy state to a coherent 
superposition of states (i.e. a "wavepacket" ) is an impor- 
tant tool in atomic and molecular physics and chemistry. 
Preparation of a wavepacket with well defined amplitudes 
and phases of the constituent eigenstates is the starting 
point for many techniques in precision spectroscopy [HH], 
coherent control of molecular dynamics and chemical re- 
actions [3], in the design of atomic clocks [I] and fault 
tolerant quantum computing schemes 5 . In many sit- 
uations, it is desirable to make the population transfer 
complete, that is to move the entire ensemble of atoms or 
molecules to the target superposition state, in a robust 
manner with respect to the excitation field parameters. 

Both adiabatic and non-adiabatic methods of complete 
population transfer between two single eigenstates are 
well established [5]. Non-adiabatic Rabi cycling in a 
two-level system enables complete population transfer by 
means of a so-called 'V-pulse " . Though a 7r-pulse can 
be realized on a very short time scale, it is far from being 
robust, as it is highly sensitive to the exact value of the 
laser intensity, frequency and pulse length. The method 
can be extended to the case of a multi-level target state 
[7J, but suffers from the following complication. 

Because of the difference in transition dipole moments 
to different excited eigenstates, a single transform- limited 
pulse can not, in general, serve as a 7r-pulse for all transi- 
tions simultaneously. Recently, we have shown that such 
non-uniformity of dipole moments can be compensated 
using the technique of spectral pulse shaping [7J. Rabi 
frequencies, and therefore pulse areas, for each transition 
can be equalized by adjusting the amplitudes of the cor- 
responding frequency components of the excitation pulse. 
The required amplitudes are, however, not easy to find 
due to the dynamic Stark shifts of the energy levels. As 
illustrated in FigjTJa) , strong polychromatic field changes 
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FIG. 1: (Color online). Illustration of the interaction of a mul- 
tilevel system with polychromatic coherent radiation consist- 
ing of a number of discrete spectral components (red vertical 
arrows). Dashed lines represent the ground state "dressed" 
with a photon. Bare excited states of the system (solid lines) 
are dynamically shifted in energy due to the presence of non- 
resonant components of the laser field. This results in in- 
complete population transfer via non-adiabatic process (a). 
When the spectral components are simultaneously chirped in 
frequency (b), adiabatic transfer into a coherent superposi- 
tion of excited states is executed. Different timing of the 
level crossings (gray arrows) complicates the control over the 
resulting wavepacket, but does not affect the completeness of 
excitation. 



the instantaneous resonant frequencies of atomic transi- 
tions, making the necessary (for a 7r-pulse ) resonance 
condition hard to satisfy throughout the whole excita- 
tion process. Similar scenario has been considered in a 
two-photon excitation scheme, where the required pulse 
shaping has been found by means of iterative feedback 
controlled loops [5]. 

In contrast to the above, the method of Adiabatic Pas- 
sage (AP) is less sensitive to the parameters of the driving 
field, transition dipole moments and resonant frequen- 
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cies [HI QUI HI]- In AP, the dynamics of a quantum 
system is described by a time-dependent Hamiltonian, 
whose instantaneous "adiabatic" eigenstates are essen- 
tially decoupled from one another as long as the field is 
strong and its parameters are changing slowly [T5J . Each 
adiabatic state is a coherent superposition of bare eigen- 
states of a field-free Hamiltonian. As the parameters of 
the interaction field change in time, projections of each 
adiabatic state on the bare-state basis set are evolving 
accordingly. By making one of the adiabatic states co- 
incide with the initial state of an atom at the beginning 
of the interaction, and with the target excited state - at 
the end, one can achieve robust and complete transfer 
of population as long as the adiabaticity conditions are 
satisfied throughout the interaction. 

Adiabatic passage between two states is most com- 
monly executed by sweeping, or "chirping" , an instan- 
taneous frequency of the driving field across a resonance. 
Adiabatic crossing of the excited state by the ground 
state dressed with a photon field results in complete pop- 
ulation transfer between those two states. AP into multi- 
level wavepackets have been recently considered in two 
distinct interaction scenarios: with multiple phase-locked 
narrowband laser fields [T3J HH US], and a single spec- 
trally shaped broadband pulse |16j . The latter approach 
provides a simple strategy of shaping the pulse in accord 
with the energy level structure of the system of inter- 
est. Introducing local frequency chirp around each res- 
onant transition frequency (Fig(ljb)), simultaneous AP 
is initiated into multiple target states which make up 
the final wavepacket. Figure [l] (b) provides an intuitive 
illustration of the robustness of the method, not only 
with respect to the field amplitude and frequency, but 
also with respect to the dynamic level shifts associated 
with the presence of a strong polychromatic field. In- 
deed, crossings of the dressed ground state with multiple 
excited states will occur even if the latter are heavily 
perturbed by the off-resonance components of the driv- 
ing field. Dynamic shifts of transition frequencies result 
in a slightly different timing of each adiabatic crossing. 
As we show here, this lack of synchronism between multi- 
ple APs slightly complicates the control of the makeup of 
the target wavepacket, but does not reduce the efficiency 
of population transfer. 

In the time domain, the interaction picture is espe- 
cially intuitive when the target states are equidistant in 
energy. In that case, shaping a pulse with multiple local 
frequency chirps results in a train of mutually coherent 
ultrashort pulses, separated by the evolution period of 
the wavepacket. Though each pulse in the train transfers 
only a small amount of population to the target super- 
position state, population of that state coherently accu- 
mulates [TJ [17j [18] piece by piece, reaching 100% at the 
end of the interaction regardless of the total energy of the 
pulse train |16| . Shown in Fig{2]is an example of numer- 
ically calculated dynamics of such Piecewise Adiabatic 
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FIG. 2: (Color online). Results of numerical simulations of 
piecewise adiabatic passage for Potassium atom. Population 
transfer is shown in (a) as a function of time. Blue line 
represents the ground state population while red and green 
lines correspond to the populations of and 4P 3 / 2 excited 

states, respectively. Excitation field amplitude is plotted in 
(b) on the same time scale. 

Passage (PAP) into a wavepacket consisting of only two 
levels and driven by a train of ultrashort pulses. 

The population transfer technique we are studying here 
can be viewed as a practical application of the 1 + TV- level 
control schemes, based on the Morris-Shore transform 
[TT?] and discussed in [T5] . The selectivity of the transfer 
is obtained by tailoring the temporal and spectral pro- 
files of the train of pulses to the target wavepacket dy- 
namics [T] [50] [5TJ [55] . The correspondence between the 
properties of the pulse train and the target wavepacket 
dynamics has also been noted in the optimization stud- 
ies aimed at either maximizing population transfer into 
a wavepacket [23] [24] [25] or stabilizing such transfer 
against the wavepacket spreading and decoherence [26] . 
Our method provides an alternative to the "multi-RAP" 
pulse sequences of [57J and "molecular 7r-pulses" used 
in exciting molecular wavepackets PS, 2S] [3U] [31] • The 
difference with the latter methods is manifest when the 
target wavepacket consists of more than two eigenstates 

OS- 

Recently, we have demonstrated experimentally the 
method of piecewise adiabatic passage into a single 
excited state [35] . Here, we extend this method to 
the multi-level target states case, by exciting the fine- 
structure doublet in atomic Potassium. We show experi- 
mentally that one can achieve complete population trans- 
fer into such "spin-orbit" wavepacket with high degree of 
robustness. We also present a simple and intuitive ap- 
proach to controlling the phases and amplitudes of the 
constituent eigenstates of the final superposition, and ad- 
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FIG. 3: (Color online), (a) Relevant quantum states of Potas- 
sium. Blue double-peak curve represents the shaped spectrum 
of pump pulses, whereas red arrows depict possible channels 
of ionizing the atoms with probe pulses, (b) Spectral shap- 
ing of pump pulses with local frequency chirping around two 
electronic resonances at 766.5 and 766.9 nm: intensity (solid 
black) and phase (dashed red). Solid green line represents 
flat spectral phase used for a non-adiabatic excitation scheme, 
(c) Experimental setup. Two Optical Parametric Amplifiers 
(OPA) are pumped by a Ti:Sapphire regenerative amplifier 
(RGA) producing 130 fs 2 mj pulses at 1 KHz repetition 
rate. One OPA is used to generate pump pulses of variable 
energy, controlled by an attenuator and measured by a fast 
photodiode. Pump pulses are shaped with a Pulse Shaper 
and are weakly focused on a cloud of Potassium atoms inside 
a vacuum chamber. Probe pulses are delayed by a variable 
time delay and tightly focused on the atomic cloud. Potas- 
sium atoms, ionized by probe pulses, are accelerated towards 
and detected by a multichannel plate based Ion Detector. 



dress its limitations. 



EXPERIMENTAL SETUP AND DETECTION 
METHOD 



The relevant levels of Potassium atom are shown in 
Figj3] (a). Two fine-structure levels, iP\/2 and 4P 3 / 2 , 
make up the target wavepacket, which is populated from 
the initial ground state AS x / 2 . The D 1 (4S , 1 / 2 — > 4P 1/2 ) 
and D 2 (4S , 1 / 2 — * 4P 3 / 2 ) transitions at 769.9 and 766.5 
nm, respectively, are excited by a single broadband laser 
pulse of 9.5 nm full width at half maximum (fwhm) cen- 
tered at 768.2 nm. The latter ("pump") pulse is pro- 
duced by a traveling wave optical parametric amplifier 
(Fig{3jc)), pumped by a Ti:Sapphire femtosecond regen- 
erative amplifier. The energy of pump pulses are at- 
tenuated and recorded with a fast photodiode prior to 
shaping. Hereafter, the reported pulse energies always 
correspond to the energy of unshaped pulses. 

To apply spectral shaping, we use a home-made pulse 
shaper implemented in 4/ geometry and based on a 



double-layer liquid-crystal spatial light modulator (SLM) 
[33] . The shaper controls both the phase and the ampli- 
tude of a pulse with spectral resolution of 0.14 nm. The 
amplitude is shaped by applying two Gaussian windows 
of variable relative height and width, each centered at 
one of the resonant frequencies (Figj3jb)) . Blocking of 
the non-resonant spectral components is important for 
achieving control over the target wavepacket, as discussed 
below. As shown in the figure, quadratic spectral phase is 
added to each Gaussian window, producing variable lin- 
ear frequency chirp around the corresponding resonance. 
Prior to applying any spectral shape described below, the 
phase of pump pulses is flattened using the technique 
of multi-photon inter-pulse interference phase scanning 
(MIIPS [31] )• Flat spectral phase is also used to drive 
Rabi oscillations as described below. The shaped pump 
pulses are then focused on a cloud of Potassium atoms 
evaporated from a Potassium dispenser inside a vacuum 
chamber. 

The population of the excited coherent superposition 
state is detected by photo-ionizing the atoms with a sec- 
ond ("probe") pulse. In order to probe the region of 
uniform pump intensity, the probe beam is focused much 
tighter than the pump pulse (beam diameters of 135 and 
460 firn, respectively). Our detection method is based 
on the weak field bichromatic coherent control scheme 
[35]. The and 4P 3 / 2 states are coupled to 5Si/ 2 

with a weak (less than 0.3 /xJ) femtosecond pulse (Fig{3] 
(a)). The central wavelength of probe pulses is tuned in 
such a way as to drive both transitions at 1243 and 1252 
nm with equal efficiency. Two more photons from the 
same probe pulse ionized the atoms. At such low probe 
energy, 5-photon ionization from the ground electronic 
state is negligible. Extracted from the interaction regions 
with a series of voltage plates, the ions are identified by 
their time of flight and detected by a micro channel plate 
detector (MCP). 

In order to understand the bichromatic control detec- 
tion we write the wavefunction of the system as: 

i 

= b {t)\AS 1/2 ) + 6i(t)|4P 1/2 ) + 6a(t)|4P 3/2 ). 

Because of the low probe energy, the ionization probabil- 
ity is proportional to the population of the 5Si/ 2 state. 
In the perturbative regime of interaction, this probability 
can be calculated as 135 : 



P(t) cx |6 1 (t)| 2 |e(o; 1 )| 2 ( i 11 + | tJ2 (t)| 2 | 
2Re[b 1 (t)b*(t)e(u; 1 )e*(u 2 )d 12 
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where &i(t) and b 2 (t) are the excited eigenstate am- 
plitudes (EqjlJ, e{ijJ\) and e(ui 2 ) are probe field ampli- 
tudes at the resonant probe transition frequencies of 
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55 1/2 and 4P 
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5 Si/ 2 j respectively, and 
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dij = {ipj\d\5Si/2){^Si/2\d\i^i) with d being the dipole 
moment operator. 

The last term in Eq.([2]) is the interference term which 
depends on the difference between the relative phase of 
the excited wavefunctions b 12 and the relative phase of 
the two resonant probe field components e(uj 1 _ 2 ). Since 
the latter is constant, the ionization probability, and 
therefore the ion signal recorded as a function of the 
pump-probe time delay, oscillates as sin[(ti>i — oj 2 )t+4>i 2 ], 
where <fii 2 = arg{6| 62^12}- These oscillations indicate 
that two 4P states are populated coherently, whereas 
their relative amplitudes and phases can be extracted 
from the oscillation contrast and phase, respectively [55] , 
In the experimental plots shown below, each data point 
is an average of 200 measurements. 

RESULTS 

We first execute population transfer into each of the 
4P states separately. Non-adiabatic interaction results 
in familiar Rabi oscillations which serve as a convenient 
tool for calibrating both the excitation pulse area and 
excited state population. The original spectrum and the 
two Gaussian windows of 1.8 nm FWHM used for am- 
plitude shaping are shown in FigQa). To drive Rabi 
oscillations on a selected transition, a single Gaussian fil- 
ter was applied around the resonant frequency of D\ or 
D 2 transition. Hereafter, we refer to the corresponding 
narrow-band pulses as "Z?i pulse" and U D 2 pulse" , re- 
spectively. In both cases, the spectral phase was kept 
flat across the whole spectrum of a pulse. Rabi oscil- 
lations are shown in Figj4] for both D x (b) and D 2 (c) 
lines. Given the available laser power, we were able to 
reach pulse areas of up to 3-7T on each transition. By fit- 
ting the data with a sin 2 A\^ 2 function, we calibrate the 
areas of D x and D 2 pulses ( A x and A 2 , respectively) ver- 
sus their energy. Simultaneously, we calibrate the mag- 
nitude of the ion signal corresponding to the pulse area 
of 7r. These calibrations are later used for assessing the 
completeness of the population transfer. 

In order to execute an adiabatic passage into a single 
excited state, a frequency chirp is introduced by apply- 
ing quadratic spectral phase shaping around one resonant 
frequency ui\. 2, i.e. tp(uj) — ^(w — lui. 2 ) 2 (Figj3^b)), while 
blocking the Gaussian window around the other reso- 
nance. The frequency chirp is gradually increased by the 
shaper until the population transfer shows AP-like satu- 
ration with intensity, such as that shown in Figsj4^b,c) by 
open circles for both D\ and D 2 transitions. Independent 
measurement of the applied chirp yields a — 270 x 10 3 
fs 2 . For both transitions, the ion signal reaches its max- 
imum at a pulse area of « ir and stays relatively flat 
with increasing pulse energy. The maximum in the AP 
efficiency is seen to coincides with a maximal number 
of Rabi oscillations, attesting to nearly complete popu- 
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FIG. 4: (Color online). Adiabatic and non-adiabatic excita- 
tion of a single (4P 1; / 2 or 4P 3 / 2 ) excited state, (a) Pump pulse 
spectra used to excite either Di or D 2 transitions separately 
(blue and red solid lines, respectively). Spectral profile of the 
initial unshaped femtosecond pulse is shown as dashed black 
curve, (b, c) Non-adiabatic (blue diamonds) and adiabatic 
(red circles) population transfer into 4P!/ 2 (b) and 4P 3 / 2 (c) 
states. Ion signal, proportional to the target state population, 
is plotted as a function of the corresponding pulse area, v4i, 2 . 
Solid lines show the anticipated sin 2 (Ai ;2 ) dependence fitted 
to the experimental data. 



lation transfer. The residual oscillations, reproduced in 
our numerical simulations (not shown), are due to the 
pixelization of the spatial light modulator. 

Once the pulse area and chirp required to satisfy adia- 
baticity conditions are determined, we can execute com- 
plete population transfer into a superposition of 4:P x / 2 
and 4P 3 /2 states. If both states are populated coherently, 
quantum beats should be observed in the ionization sig- 
nal as a function of the probe pulse delay according to the 
interference term in Eq.Q. In both the non-adiabatic 
and adiabatic approaches, the amplitude of the spectrum 
of the pump pulses is shaped by opening both the D\ and 
D 2 Gaussian windows simultaneously. Interference of the 
D\ and D 2 pulses in the time domain results in a train 
of pulses separated by a period of the quantum evolution 
of the wavepacket. The latter is inversely proportional 
to the fine-structure splitting of 4P state (1.73 THz) and 
equals 578 fs. A numerical example of such pulse train is 
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FIG. 5: (Color online). Non-adiabatic population transfer 
into a superposition of 4Pi/2 and 4P3/2 states. The two- 
dimensional plots (a: experiment, c: calculation) show the ion 
signal (color coded) as a function of the effective pulse area 
and pump-probe time delay. Oscillations along the vertical 
(time) axis reflect quantum beating due to the time evolution 
of the wavepacket. The quantum beats are shown in (b) for 
pulse areas of tt (blue) and 2n (red), corresponding to the 
dashed vertical white lines in (a). In plot (d), we present the 
calculated populations of the two excited states, ^P\/2 (blue) 
and 4P 3 /2 (green), and the relative phase (red) between the 
corresponding wavefunctions, as a function of the effective 
pulse area of the excitation field. The ground state population 
is shown as a dashed black line. 



shown in Fig{2^b). 

First, the spectral phase of pump pulses is kept flat 
across the whole double-peak spectrum. In the time do- 
main, flat spectral phase translates into a train of pulses 
with constant carrier oscillation frequency and no ex- 
tra phase shift between consecutive pulses [16]. The re- 
sulting Rabi oscillations between the ground state and 
the excited wavepacket are shown in Fig[5|a). The two- 
dimensional plot shows the ion signal as a function of 
the time delay (vertical axis) and effective pulse area 
(horizontal axis), which in the case of multiple excita- 
tion channels can be used as a convenient scale of the 
interaction strength. It is defined as: 



A eS = JA\ + A\ 



(3) 



where Ai = fli(t)dt, with Qi being the time depen- 
dent Rabi frequencies for the i-th transition. 

Oscillations along the vertical axis indicate quantum 
beating between 4P!/ 2 and 4P 3 / 2 states. Both the ex- 
perimental and numerical plots (Fig{5^a,c)) demonstrate 
strong dependence of the result of the excitation on the 
effective pulse area. Rabi oscillations between the ground 
state and the excited wavepacket are manifested by the 



periodic re-appearance of the beat signal. Thus, high 
beating contrast at A h = tt disappears almost com- 
pletely at A c g ss 2ir, as shown by the vertical cross- 
sections of the two-dimensional data (white dashed lines) 
in plot (b). 

We note that both the contrast of the beat signal 
and its phase are affected by an increase in the pulse 
energy. We attribute this phase change to the energy- 
dependent AC Stark shifts of the levels, caused by the 
strong non-resonant components of the polychromatic ex- 
citation field (see schematic illustration in Fig {I]). In the 
case considered here, the 4P 1 / 2 level is shifted due to the 
presence of an off-resonant Z3 2 pulse, whereas the 4P 3 / 2 
level is shifted by an off-resonant D\ field. Since the shifts 
are in general unequal, the accumulated quantum phase 
of the two wavefunctions, |4Si/ 2 ) and |4S 3 / 2 ), depends 
on the energy of both pulses. The effect is reproduced 
by the numerical calculations shown in Figj5^c). 

For pulse areas of order tt and higher, typically needed 
for significant population transfers, such dynamic cross- 
talks between simultaneously driven transitions become 
substantial when the pulse bandwidth approaches the en- 
ergy separation between the levels. In this case, dynamic 
Stark shifts cannot be neglected and multiple interac- 
tion channels cannot be treated independently. This sig- 
nificantly complicates the non-adiabatic dynamics, mak- 
ing controlled and complete population transfer hard to 
achieve. Figjf^d) demonstrates the degree to which the 
dependence of the |&i(£)| 2 and |^2(i)| 2 populations on 
pulse energy deviates from the periodic Rabi oscillations 
behavior, once the effective pulse area exceeds tt. 

In striking contrast to Rabi flopping, when local chirp 
is added to both the D\ and D 2 spectral windows (Fig. 
|3jb)), the observed quantum beats become insensitive 
to the pulse area (Fig|6ja,c)). Using our experimental 
observation of separate adiabatic passages into each level 
(Figj4|, we set the frequency chirp to 270 x 10 3 fs 2 , for 
both the Di and D 2 pulses. As seen in Figj6^b), the 
contrast of the experimentally observed quantum beating 
changes only little throughout the wide range of pulse 
areas, A e g w tt to A c ff rs 3tt. This demonstrates the 
stability of the population transfer against pulse energy 
in agreement with the adiabatic passage scenario. 

Though clearly more stable than in the non-adiabatic 
regime, the beat signal gives only indirect evidence of the 
measure of the robustness of the population transfer. Un- 
fortunately, extracting the absolute values of the popu- 
lation transferred to 4P\/ 2 and 4P 3 / 2 from the measured 
beat signal proved inaccurate. We therefore calculate 
these populations numerically. The results, presented in 
Figj6^d), show that the population ratio varies with the 
excitation pulse energy. A number of simple approaches 
to controlling this ratio are discussed later in the text. 
For the time being we note that an increase in the pulse 
area leaves the phase of the excited wavepacket intact. 
This statement is confirmed experimentally and theoret- 
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FIG. 6: (Color online). Adiabatic population transfer into a 
superposition of 4Pl/2 and 4P 3 / 2 states. The two-dimensional 
plots (a: experiment, c: calculation) show the ion signal 
(color coded) as a function of the effective pulse area and 
pump-probe time delay. The oscillations along the vertical 
(time) axis reflect quantum beating due to the time evolution 
of the wavepacket. The quantum beats are shown in (b) for 
pulse areas of tv (blue) and 2n (red), corresponding to the 
dashed vertical white lines in (a). In plot (d), we present 
the calculated populations of the two excited states, 4P 1( /2 
(blue) and 4P 3/ / 2 (green), and the relative phase (red) be- 
tween the corresponding wavefunctions, as a function of the 
effective pulse area of the excitation field. The ground state 
population is shown as a dashed black line. 

ically, as shown in panels (b) and (d), respectively. 

In contrast to the non-adiabatic dynamics, AP exhibits 
a threshold at A c g ss it, beyond which the ground level 
remains empty, even though the population ratio be- 
tween the excited levels may vary. This stability of the 
completeness of population transfer against variations in 
pulse energy is best demonstrated by comparing the nu- 
merical results of Figs J5]^d) and[6jd). It also forms the 
basis for the proposed schemes of controlling the popula- 
tion distribution among excited states. We therefore seek 
an independent experimental confirmation of the com- 
pleteness of population transfer. 

In the experiments described above, we have used Rabi 
oscillations between the ground state and one of 4P ex- 
cited states for calibrating the population transfer effi- 
ciency, assuming that the first maximum of these oscilla- 
tions corresponds to a complete transfer. However, even 
though a 7r-pulse is expected to drive the whole popula- 
tion to the excited state, it is conceivable that complete- 
ness is not fulfilled due to the sensitivity of the latter to 
frequency detuning. In order to perform an independent 
check of the completeness of the population transfer, we 
have carried out the following experiment: We delay the 
D\ pulse by approximately 8 ps with respect to the D2 
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FIG. 7: (Color online). Demonstration of the population 
transfer completeness, (a) Ion signal as a function of the 
pump-probe time delay for a sequence of frequency chirped 
D\ and D2 pulses separated by 8 ps (marked by dashed ver- 
tical line). Blue (red) dots correspond to the signal before 
(after) the arrival of the second pulse, (b) Fourier transform 
of the ion signal before (blue) and after (red) the arrival of 
the second pulse (log scale). Vertical dashed line marks the 
frequency of the quantum beating (1.73 THz) where a strong 
peak is expected in the case when both states are populated. 

pulse by inserting a thin piece of glass in front of the cor- 
responding spectral window in the Fourier plane of the 
pulse shaper. If the D2 pulse, arriving first, drives all 
the atoms to P 3 / 2 state and depletes the ground state, 
the D\ pulse, arriving second, will leave the system un- 
changed. As all atoms are residing in a single P3/2 state, 
no quantum beats would be observed. Alternatively, if 
the excitation driven by the first pulse is not complete, 
the second pump pulse would move the remaining ground 
state population to P1/2 state, resulting in the appear- 
ance of quantum beating. 

To verify this, we set the areas of both D\ and D2 
pulses (now separated in time) to 7r, and their frequency 
chirp to 270 x 10 3 fs 2 . As can be seen in Figsj4^b,c), these 
parameters are sufficient for reaching maximum popula- 
tion transfer to either P1/2 or P3/2 when the correspond- 
ing pulse is acting alone. The results of the interaction 
with a sequence of two pulses are presented in Fig|7] The 
first pulse arrives at approximately ps, while the time 
of arrival of the second pulse is about 8 ps. Fourier anal- 
ysis of the ion signal after the arrival of the second pulse 
(Fig{7jc)) shows essentially no quantum beating, ensur- 
ing that the population of 4P 1( /2 state is below 3%. 

As pointed out by us in the past [16], piecewise adi- 
abatic passage potentially combines the efficiency and 
robustness of AP with the ability to excite complex su- 
perposition states (e.g. wavepackets) and control the 
makeup of the excited state wavefunction [IB]. Control 
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over the phase of a spin-orbit wavepacket in Potassium 
has been previously demonstrated experimentally in the 
weak field regime [37]. Here we show that the relative 
phase between the two eigenstates of the target super- 
position can be controlled even when the field is strong 
enough to ensure the adiabaticity of process. Control 
over the phase is implemented by adding an extra con- 
stant phase to one of the pulses, D\ or D 2 , on top of 
the local frequency chirp in the respective spectral win- 
dow. The resulting phase profile, attained using the 
pulse shaper, is shown in Fig{8^a) for a relative phase 
of 0, 7r/2 and tt radian. Two-dimensional energy-time 
scans for these three phase shifts are shown in panel (b). 
The vertical cross-sections of each two-dimensional plot, 
displaying quantum beat patterns corresponding to » tt 
pulses, are plotted in panel (c). The phase of the oscil- 
lations, which directly reflects the relative phase of P\ji 
and P3/2 eigenstates, closely follows the extra phase shift 
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FIG. 8: (Color online). Experimental control over the quan- 
tum phase of an excited wavepacket. (a) To control the rela- 
tive phase of the eigenstates in the target superposition state, 
an extra constant phase shift of (red), 7r/2 (green) and n ra- 
dian (blue) is added to the spectral phase of D\ pulses. Black 
dashed line shows the double-peaked spectral amplitude of 
the pulse. The introduced phase shift results in a correspond- 
ing vertical shift of the detected quantum beats shown in (b) . 
(c) Vertical cross-sections of the two-dimensional scans in (b) 
at pulse area of 7r. 
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FIG. 9: (Color online). Quality of population control with 
piecewise adiabatic passage. Numerical simulations demon- 
strate possible discrepancy between the achieved population 
ratio (3 of two AP states (vertical axis) and its target value 
(horizontal axis), for different control methods (see text for 
symbol description). 



introduced via the pulse shaper. 

One may expect that the amplitudes of the eigenstates 
in the target superposition can be controlled individu- 
ally, similarly to the phase control described above. As 
was shown in [16] for the case of negligible Stark shifts, 
the state amplitudes at the end of the interaction are sim- 
ply proportional to Rabi frequencies of the corresponding 
transitions. Thus, the distribution of populations among 
excited states (here, 4Pin and 4P 3 / 2 ) can be controlled 
by changing the relative strength of the corresponding 
spectral components of the excitation field (here, the en- 
ergies of D\ and D 2 pulses). This simple strategy fails, 
however, once the dynamic Stark shifts become compa- 
rable to the energy bandwidth of the laser pulses. This 
situation is illustrated in FigjIJb), where one can see that 
unequal Stark shifts effectively change the time at which 
each adiabatic passage is executed. The farther these 
APs from being simultaneous, the bigger the deviation 
of the population distribution from the expected one. 

Poor accuracy in retrieving state populations from the 
observed beat signal does not allow us to investigate the 
possibilities of amplitude control experimentally. Here, 
we present numerical analysis of various ways of achieving 
reasonably high degree of control over the excited state 
populations by either avoiding, or compensating for the 
detrimental effects of the Stark shifts. The results are 
summarized in Fig]9] in which the calculated ratio (3 = 
\bi\ 2 1 1 62 1 2 is plotted versus its target value for different 
control methods. In all calculations, the effective pulse 
area was set to w, and the local spectral chirps near D\ 
and D 2 lines was equal to 270 x 10 3 fs 2 . 
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We first note that for the parameters used in our ex- 
periments, calculated values of p (red circles) lie within 
37.5% of those expected from a simple model which does 
not take Stark shifts into account (diagonal dashed line). 
Assuming the correct phases of the state amplitudes, this 
corresponds to the transfer fidelity (projection on the tar- 
get wavepacket) varying between 1 and 0.997. The dis- 
crepancy can be further reduced by narrowing the spec- 
tral bandwidth of the excitation pulses. Decreasing it 
from 1.8 nm (as used here) to 0.18 nm, dramatically im- 
proves the quality of amplitude control, marked as blue 
diamonds. It is important to note that such improve- 
ment comes at the expense of much longer (by a factor 
of 10) interaction time, and higher spectral resolution of 
a pulse shaper, which will be necessary for the frequency 
chirping of a spectrally narrower pulse. 

Unlike the case of independent adiabatic passages, i.e. 
when Stark shifts can be ignored, the Stark-induced dy- 
namical cross-talk between different AP channels makes 
(3 sensitive not only to the ratio of D\ and D2 pulse en- 
ergies, but also to their sum. This certainly reduces the 
robustness of the proposed method of amplitude con- 
trol, though not the robustness of PAP itself. In other 
words, in PAP, the population transfer remains complete, 
even though the shape of the excited wavepacket may 
change in response to changes in the overall pulse energy. 
Utilizing this important property of piecewise adiabatic 
passage, we suggest a simple adaptive strategy of control- 
ling the eigenstate amplitudes in the excited superposi- 
tion state. By accurately measuring populations of the 
excited states (not available in our current experimental 
setup), a correction can be introduced into the relative 
energy of each spectral component of the excitation field, 
proportional to the deviation of the observed population 
distribution from its target shape. The numerical results 
of such an iterative procedure appear as green triangles 
in Fig(9] After only two iterations, the calculated ratio 
between |&i| 2 and | 1 2 is brought within 0.2% of its tar- 
get value! Unlike the more general schemes of adaptive 
control [38 , the present technique involves only as many 
control variables as there are excited states (in our case, 
two), and converges very quickly due to the inherent ro- 
bustness of AP. 



CONCLUSION 

We have presented experimental and numerical stud- 
ies of the population transfer using piecewise adiabatic 
passage from a single ground state to a superposition 
of excited states. As in conventional adiabatic popula- 
tion transfer into a single state, our method of piecewise 
adiabatic passage into a wavepacket is insensitive to the 
driving field amplitude. The method allows for complete 
population transfer and offers control of both the phase 
and amplitudes of the eigenstates composing the target 



superposition state. The latter are retrieved by applying 
bichromatic control to the process of photoionization. 

This work has been supported by the CFI, BCKDF, 
NSERC and DTRA. 
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